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Abstract 


A hybrid method that combines the finite element method (FEM) 
and the boundary element method (BEM) is developed to analyze elec- 
tromagnetic scattering from arbitrarily shaped material cylinders. By 
this method, the material cylinder is first enclosed by a fictitious 
boundary. Maxwell's equations are then solved by FEM inside and by 
BEM outside the boundary. Electromagnetic scattering from several 
arbitrarily shaped material cylinders is computed and compared with 
results obtained by other numerical techniques. 

1. Introduction 

The problem of electromagnetic scattering determination from an arbitrarily shaped material cylin- 
der is considered in this paper. Electromagnetic scattering from an inhomogeneous material cylinder 
can be determined by two basic approaches — an integral equation and a differential equation. 

In the integral equation approach, a surface integral equation or a volume integral equation may be 
used. In the surface integral equation formulation, conducting surfaces of a scatterer are replaced by 
equivalent surface electric currents, whereas material surfaces are replaced by both equivalent surface 
electric and magnetic currents. Coupled integral equations are then formed by application of appropriate 
boundary conditions to the different field components that are produced by the equivalent currents. 
These coupled integral equations are solved by the method of moments for the unknown equivalent cur- 
rents. Arvas and Sarkar (ref. 1) have analyzed radar cross sections of various two-dimensional struc- 
tures with the surface integral equation approach. 

For a cylinder consisting of inhomogeneous material, the surface integral equation approach does 
not correctly model the inhomogeneity of the material cylinder. In such cases the volume integral equa- 
tion formulation is then used to accurately model the inhomogeneous material cylinder. In the volume 
integral equation formulation, the material cylinder is divided into unit cells. These unit cells may be 
rectangular bricks or tetrahedral shapes and are small enough so that the field intensity is nearly uniform 
within each cell. The material cylinder is then assumed to be replaced by an equivalent polarization cur- 
rent flowing in these unit cells. The coupled integral equations, which are formed by application of 
appropriate boundary conditions to the field produced from the equivalent currents, are then solved by 
the method of moments. Using the polarization current concept, Richmond (ref. 2) determined scattered 
field patterns of dielectric cylinders of arbitrary cross-sectional shape. The equivalent current concept 
has also been used by Sarkar and Arvas (ref. 3) and by Schaubert, Wilton, and Glisson (ref. 4) to ana- 
lyze electromagnetic scattering by arbitrarily shaped inhomogeneous dielectric bodies. However, the 
method involves the solution of a fully dense matrix equation, which may require prohibitively large 
computer memory and long computation time. 

In the differential equation approach, the volume of the inhomogeneous scatterer is discretized, and 
Maxwell’s equations in differential equation form are solved. Some advantages of the differential equa- 
tion approach over the integral equation approach are: (1) a method more convenient for handling com- 
plex inhomogeneous scatterers, (2) a method that results in sparse matrix equations, and (3) a method 
better suited for closed-region problems. However, for open-region problems in electromagnetics, such 
as radiation and scattering, the volume must be properly tmneated with an artificial or fictitious bound- 
ary with proper absorbing boundary conditions. The accuracy of results obtained with an artificial 
boundary with absorbing boundary conditions depends upon boundary locations as well as the order of 
the boundary conditions. 

An alternative approach (known as the hybrid approach) presented in this paper retains the advan- 
tages of both differential equation and integral equation approaches. The general procedure for a hybrid 
technique requires that the scatterer be enclosed by an artificial boundary. Maxwell’s equations are then 



solved by a differential equation approach such as the finite element method (FEM) inside the artificial 
boundary and by an integral equation approach in discretized form such as the boundary element 
method (BEM) outside the artificial boundary. Although use of BEM on and outside the artificial 
boundary results in a full dense matrix, convergence of an approximate solution to the exact solution is 
guaranteed without a change in the location of the artificial boundary. 

The remainder of the paper is organized as follows. Section 3 contains basic formulation of the 
problem for FEM and BEM. Scattering amplitudes of a few sample problems are computed by the 
present technique in section 4. Also in section 4, results obtained by the present method are compared 
with those obtained by other analytical techniques. The paper is concluded in section 5 with comments 
on the future scope and extension of the present work. 

2. Symbols 

A € jj (I, J)th element of FEM matrix for eth subdomain 

a cylinder radius; triangle side length 

[B] finite element matrix 

b boundary circle radius; triangle base length 

C fictitious boundary curve 

[. D ] sparse matrix 

D Jp (J, p)ih element of matrix [D] 

Eq amplitude of incident electric field 

E z Z- axis component of total electric field 

Em incident electric field vector 

e Jwt time convention 

/ frequency, Hz 

f s E or H scattered far field 

G(.) free-space Green’s function, where (.) represents appropriate argument 

[G] square matrix of order P xP 

Gpq (p, q ) th element of matrix [G] 

[G] _1 inverse of matrix [G] 

[H] square matrix of order P X P 

H pq (p, q) th element of matrix [H] 

( 2 ) 

H l ( . ) Hankel function of first order and second kind 
( 2 ) 

Hq (.) Hankel function of zero order and second kind 
H z Z - axis component of total magnetic field 

Hq amplitude of incident magnetic field 

H in incident magnetic field vector 

/, J integers 

f z equivalent electric current representing incident field 

j imaginary number equal to V— 1 

propagation constant in free space 
m, n integers 
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p 

p 

Q p (r 1 ) 
q 

r 

T{x , y) 

U P 

w] 

( x 2 ’ ^ 2 )’ ' 

( x 3’ ^ 3 ) , 

(*> y) 
z 

a r (x,y), 
$ r (x, y) 

[H 

[H c 

r. 


* P - 1 
*// 

e(x, y) 
£ r (x, y) 


ri 

e P 

H(x,y) 
P r (x, y) 
V 0 

(p,<l>) 
(p', f) 

G 2D 

m c 

Y « 

*t*in 


unit normal vector of curve C drawn outward 
unit normal vector of curve C drawn inward 
number of nodal points on curve C 
integer for pth nodal point 
piecewise distribution along r| direction 
integer for ^th nodal point 
unit vector 

testing function equivalent to W e f 

column matrix with nodal amplitudes of normal derivative on curve C 

amplitude of normal derivative of \\f t on curve C 

expansion function for 7th node over eth triangular subdomain 

coordinates of three vertices of triangle 

variables of rectangular coordinates 
unit vector along Z-axis 

equal to p r (x, y), 8 r (x, y), respectively, for TM case; equal to e r (x, y), p r (x, y), 
respectively, for TE case 

column matrix for nodal amplitude 
column matrix for nodal amplitude on curve C 
amplitude of at 7th node 

length of linear segment between p and p + 1 nodes on curve C 
length of linear segment between p and p — 1 nodes on curve C 
Kronecker delta function 

e 0 M x ’ y) 

relative permittivity at (x, y) 

permittivity of free space 

variable along linear segments of curve C 

angle between two consecutive linear segments of curve C, deg 

PoM x ’ y) 

relative permeability at (x, y) 

permeability of free space 

variables of cylindrical coordinate system 

source coordinate 

radar scattering width 

column matrix with nodal amplitudes of incident field on curve C 
amplitude of incident field at qth node on curve C 
incident angle of electromagnetic wave 
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\|/ in incident electric or magnetic field scalar 

\\r t total electric field E z for TM case; total magnetic field for TE case 

co angular frequency ( 2%f ), rad/sec 

Abbreviations: 

in incident 

TE transverse electric 

TM transverse magnetic 

3. Theory 

Consider a material cylinder infinite along the Z-axis and of arbitrary cross section as shown in fig- 
ure 1 . For the purpose of analysis, the scattering structure is divided into regions I and II by contour C. 
Region II may consist in general of inhomogeneous material with permeability |LL(x, y) and permittivity 
e(x, y). It may also have embedded metallic strips within the material. Region I is free space with perme- 
ability p,() and permittivity Eq surrounding the cylinder. An incident electromagnetic wave is assumed in 
the direction normal to the axis of the cylinder. Both transverse electric (TE) and transverse magnetic 
(TM) polarizations are considered, and e JG>t time convention is assumed. 

For TM and TE incidence, respectively, the incident field is given by 


yv j-, ^oPCOS^-^in) ] 

in = z¥ . n = z E q l 

J k oP cos ((j) — <j> in ) [ 

in — ZVin — 


(1) 


where E 0 = H 0 = 1 is the amplitude, <{) in is the incident angle, and (p, (J)) are the variables of the 
cylindrical coordinate system. The total field \) f t inside the region bounded by curve C can be deter- 
mined from the solution of the source-free scalar wave equation given by 


V- 


La r (*, y) 




+ * 0 P r (*» = 0 


( 2 ) 


In equation (2) a r (x, y ) , [3 r (je, y ) , and \|/ f are equal to p r (jc, y) , £ r (x, y) , and E , respectively, for TM 
Z-axis excitation and are equal to £ r (x, y) , p. r (x, y) , and H , respectively, for TE to Z-axis excitation. 
Multiplication of both sides of equation (2) by a testing function T(x, y) and integration over the cross 
section of region II yields a weak form of the wave equation (Silvester and Ferrari, ref. 5) as 


J J r | V ‘ [ a (xTy) 7 ^* + P'-M*’ [ dx dy = 0 

Region II ^ r * 


With the use of vector identity 

1 




la r (x, y) 




= TV • V 


1 


a r (x, y) 


Vt 


+ vr- v 


l 


_a r (x, y) 


% 


(3) 


(4) 


and the divergence theorem, equation (3) can be written as 


I i 


Region II 


vr - v 


l 


la(x, y)^J 


- & 0 P r (.x, y)T\ \f { > dx dy = ^rV\j/ ? • h dc 


where h is a unit vector drawn outwardly to the curve C as shown in figure 1. 


( 5 ) 
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To construct an approximate solution of equation (5) by FEM, region II is approximated by a union 

of triangles as shown in figure 2(a). On the eth triangle \|r is represented by a linear combination of 
€ * 
functions Wj(x,y) as 


where 


3 

\| f e { = ^TjW^x, y) (on eth triangle) 
1=1 

\j/^ = 0 (otherwise) 


W]{x,y) = 


1 X 1 y x 

-1 

1 ~ 

1 x 2 yi 


^12 

1 *3 ?3_ 


j>/3 


( 6 ) 


and the Kronecker delta function is defined as 5 u = 0 for 7^7 and = 1 for 7 = 7. The ampli- 
tudes Tj, Y 2 , and r 3 are the unknown values of \j/ t at the three vertices of the eth triangle having coordi- 
nates (xj, y x ), (x 2 , y 2 ), and(x 3 ,y 3 ), respectively. By Galerkin’s method with T(x, y) = Wj(x, y) and 
7=1,2, and 3, the left side of equation (5) over the eth triangle can be written as 




1=1 


VW r ■ V 


L« r (x, y) 


W, 


-k$ r {x,y)WjWj \dxdy (7=1,2, 3) 


( 7 ) 


The previous expression can be written in a matrix form as 


where 



l u 


II 


VW • V 


2 w e 

La r (x, y) 1 


k 2 t fi r (x,y)W e J W, e 


dx dy 


( 8 ) 


( 9 ) 


Now consider a union of two triangles as shown in figure 2(b); an expression corresponding to 
equation (8) for the union of two triangles can be written as 


A 6 

A 11 

A e 

a 12 

A e 
A 13 

1 

O 


r l 

A e 

A 21 

(4 2 + 4s) 

(4* + 4>) 

A f 


r 2 

A e 

a 31 

(4l + A f 65 ) 

(A 33 + Ag 6 ) 

A f 

a 64 


r 3 

1 

O 

. j 

a 45 

A f 

n 46 

0 


r 4 


f 

where A u are determined from equation (9) with superscript e replaced by superscript /and the integra- 
tion performed over the / triangle. An appropriate alignment of common edges, as indicated in 
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figure 2(b), is required to assure the correct combination of terms in the matrix equation (10). To ensure 
continuity of fields across the common edge, T 6 = T 3 and T 5 = T 2 were enforced in the derivation of 
equation (10). With an assembled mesh of all triangles over the surface bounded by curve C, the left 
side of equation (5) can be written as 


[B\[ T] (11) 

where the elements of the matrix [B] are obtained from [T] is the column matrix with elements 

u J Win 

given by the values of \|/ f at the vertices of triangular elements, m — 1, 2, 3, ..., N, n = 1, 2, 3, ..., N, and 
N is the total number of nodes. 

To evaluate the contour integration on the right side of equation (5), the normal derivative of \j/ r 
over curve C is required. To determine the normal derivative of \|/ f on curve C, the following procedure 
is used. 

The boundary curve C is discretized into linear segments as shown in figure 3 and the normal deriv- 
ative of the function \j i t on the curve is written as 


P 

V V ,« = ££/ p e p (Tl) (12) 

p = 1 

where U p is the unknown amplitude at pth node on curve C. The function Q p (j\) varies linearly with T| 
over the segment as 


e p (iD = t^ 1 - (o< t| < a ,) 

p-1 

A -(Tl-A l} (13) 

Q P W=- A P (Ap_ i — T1 - A^) 

P J 

The right side of equation (5) can therefore be written as 

p 

§TVv/ t hdc = ^U p j T Q p m)dr\ = [D][m (14) 

C p=l c 

where the elements of matrix [£)] are given by 

d j p = §KQ P ^ di \ ( 15 > 

C 

The unknown amplitude U p is determined as follows. The function \|/ ? in region I is obtained by the 
solution of the inhomogeneous wave equation 

V 2 y t + k 2 0 \ \f t = -j(x>VL 0 J l z (16) 

subject to essential and natural boundary conditions on curve C. These essential and natural boundary 
conditions are the values of the function and its normal derivative on curve C. In equation (16), f z is the 
equivalent current source producing the incident field. 
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To obtain a solution, multiply equation (16) by G(p, <|>; p', (j)'), where G(.) is the free-space Green’s 
function for a line source located at (p', 6'), and integrate over region I, which gives 


J | ( V2y, + $|f,)G(.) dx dy = -7'0J(X 0 j J f z C(.) dx dy 

Region I Region I 

From the scalar Green’s theorem 

J | [G(.)V 2 \\f t -\\f t V 2 G(.)] dx dy = hi dc ~ j>\\f { VG(.) • h { dc 

Region I C C 

where the unit normal vector h j is shown in figure 1. Equation (17) can then be written as 
J lv,lV 2 G(.) + k 2 0 G(..)]dxdy = -jw^ J jj‘ z G(.) dx dy 

Region I Region I 


+ o\|/ r VG • dc — <jJV\|f, • n x G dc 


(17) 


(18) 


(19) 


Because [V 2 G(.) + &qG(.)1 = — 6(p — p')5(<|) — <j)')/p, j'cop J JV^G(.) dx dy = \|/ f , 

Region I 

= -h, equation (19) can be written as 

i|/ ? (p, (()) = V,-(p, <|>) - <j>V\|/ r • fiG(.)d'X) + <J)\|/ r VG(.) • h di] 
c c 

The free-space Green’s function used in the previous equations is given by 

G (■) = jj.H^olp-pi) 

VG(.) = ^H ( , 2) (* 0 IMl) 


and 


( 20 ) 


(21) 


( 22 ) 


( 2 ) ( 2 ) 

where H 0 (.) and H\ (.) are the Hankel functions of the second kind of zero and first order, respec- 
tively, and r is a unit vector along p - p'. The function \| f t on the curve in figure 3 can now be repre- 
sented by 


v, = p 2p(n) 

P = 1 


(23) 


where T p is the unknown amplitude of \\f t at the pth node on curve C. Substitution of equations (12) 
and (21)-(23) into equation (20) yields 

'k P 


p - 1 c 


{Yd U p\ H o\o\h-9'\)Q p Wdy\ 
p - 1 c 


(24) 
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where is the value of incident field at the qth node on curve C. Equation (24) can be written in 
matrix form as 


where 


and 


H n Hn ■ 
H 2\ H 22 • 

. H lp 

. h 2P 


r 2 


*1 

— 

G n g \2 • 
G 2\ G 22 • 

■ G \P 

. g 2P 

H P\ H P2 • 

■ H PP 


?e 


Jp 


G P\ G P2 • 

• G pp 


u 2 

Vp_ 


H. 


pq 




qq 


(25) 


(26) 


G Pq = ijH^(k 2 0 \ Pll -p’\h p (r)) dt\ (27) 

c 

The normal derivative of \|/ ? is then obtained by solution of equation (25) for [U] c in terms of [r] c to 
get 


[ u] c = [G] ‘[//Kric-IG] '[T] c ( 28 ) 

where [U] c and [r] c refer to unknown nodal amplitudes on curve C. With substitution of equa- 
tion (28) into equation (14), the rightside of equation (5) reduces to 

|rv«i»,«d C = [ojiGr'itfHnc-iDHGr'm,.. (29) 

c 

From equations (1 1) and (29), equation (5) can be written in matrix form as 

m[n-[£>][GrV][r] c = [DKGr'm,, (30) 

and solved for [T], which also includes additional [T] c terms. Then [ U] c is determined by substitu- 
tion of [r] c into equation (28). The scattered far field in the cp direction is then determined by the last 
two terms of equation (24) with = p and the asymptotic evaluations of the Hankel functions as 


, k 0( 2 j V /2 ~ jk oPVr ^ Jk o P'cos(9-(p') 

^ = e X r Je p (n) e <t*\ 

p = 1 C 


Jt&oP, 


(31) 


p = l c 
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where (p', cp') is a function of r|. For TM and TE incidence, respectively, the radar cross section is then 
obtained from 


°2 D 
a 2 D 


limit p oo 2 tc p 

limitp^ 2k p 




(32) 


4. Numerical Results 

In this section, numerical results for the scattering width of various 2-D material cylinders are pre- 
sented. To solve equation (30), the matrices [5], [D], [G], and [ H ] must be determined. Matrices 
[2?] and [£>] are sparse matrices, and their elements are determined by equations (9) and (15), respec- 
tively. Matrices [H] and [ G ] are dense matrices, and their elements are determined by equations (26) 
and (27), respectively. Equations (26) and (27) are evaluated with the use of Gauss quadrature numeri- 
cal integration. When p = q, the integration in equation (26) results in H qq = (Q q /(2n)), where 0^ is 
the internal angle at the qth. node shown in figure 3. The column matrix [F] , which is obtained after the 
solution of equation (30), is used in equation (28) to determine [U] c . From known [F] and [t/] c 
on curve C, the scattered far field, and hence, the scattering width are determined by equations (31) 
and (32). 

To validate the computer code, the bistatic scattering width of a perfectly conducting cylinder of 
radius 1 .OX, when excited by a TM-polarized plane wave, is calculated by the present formulation and 
compared with the calculated results given in reference 1. (See fig. 4.) Results of the two calculation 
methods are in good agreement. For calculations by the present method, the circular cylinder is assumed 
to be enclosed by an artificial circle with radius equal to 1 .2X. The region enclosed by the artificial cir- 
cle is analyzed by the finite element method, and the region outside the circle is analyzed by the bound- 
ary element method. For the BEM the circle is divided into 76 points. Selection of the location of the 
artificial boundary around the cylinder is arbitrary. For the conducting cylinder, when the location of the 
artificial boundary coincides with the object boundary, the problem can be solved by the use of only 
BEM. For validation of the FEM formulation, the location of the artificial boundary was selected at 1.2X 
radius to provide a minimum one- cell-thick FEM region. If the artificial boundary is moved farther 
away from the object boundary, the computational area is unnecessarily increased. The circumference 
of the artificial boundary was divided into 76 linear segments so that the length of each small segment 
was less than 0.1X. In figure 5 the bistatic scattering width of the circular cylinder, when excited by a 
TE-polarized plane wave, is shown and compared with earlier published results. (See ref. 1.) 

To check if dielectric and magnetic materials are properly handled by the present method, the 
bistatic scattering width of a material-coated circular cylinder is calculated for both TE and TM excita- 
tions and presented in figures 6-8 along with earlier published results. (See ref. 1.) Overall earlier 
results and the present method results are in good agreement. For the results presented in figures 6-8, 
the artificial circular boundary was taken at radius b. The artificial boundary was divided into 90 seg- 
ments. Any further increase in the discretization of the artificial boundary C resulted in insignificant 
changes in the results presented in figures 6-8. 

In general, the artificial boundary curve C enclosing the cylindrical scattering structure may be of 
any shape. To check this aspect, the bistatic scattering widths of cylinders with triangular cross sections, 
as shown in figure 9(a), are considered. The triangular cylinders are enclosed by the following bound- 
aries: (1) circular boundary (fig. 9(a)), (2) elliptical boundary (fig. 9(b)), (3) conformal boundary with 
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blended comers (fig. 9(c)), and (4) conformal boundary with sharp comers (fig. 9(d)). The bistatic scat- 
tering widths, which were obtained for the four examples, are shown in figure 10 for comparison with 
earlier results published by Peterson and Castillo in reference 6. Results presented in figure 10 are in 
good agreement with each other. The results presented in figure 9 show advantages from the artificial 
boundary selected to be as close as possible to the outer boundary of the cylinder. This reduces the com- 
putational surface and, hence, the number of unknowns. For the examples considered, the artificial 
boundary selected in figure 9(d) is the optimum because it has the least number of unknowns when 
compared with the artificial boundaries considered in figures 9(a)-9(c). 

Other geometries and their results are shown in figures 11-13. The geometries that were considered 
included the conducting strip shown in figure 11(a), the microstrip transmission line shown in 
figure 12(a), and the conducting cylinder of von Karman shape shown in figure 13(a). For the calcula- 
tion of the scattering widths of a conducting strip, a microstrip transmission line, and a conducting cyl- 
inder of von Karman shape, the artificial boundaries were selected to be elliptical enclosures of the outer 
boundaries of the structures as shown in figures 11(a), 12(a), and 13(a), respectively. The scattering 
widths of these structures calculated by the present method are presented along with earlier published 
results (ref. 1) in figures 11(b), 11(c), 12(b), and 13(b). The results of all geometries are in good 
agreement. 

5. Conclusion 

A hybrid technique that combines the finite element method and boundary element method has been 
used to determine backscattered fields from arbitrarily shaped material cylinders. Validity of the com- 
puter code developed with the present hybrid technique has been demonstrated with various arbitrarily 
shaped material cylinders. Although use of the boundary element method in the present hybrid tech- 
nique leads to a partly sparse and partly dense matrix, the present method is guaranteed to converge irre- 
spective of the shape of the terminating or artificial boundary enclosing the arbitrarily shaped material 
cylinder. The present hybrid technique can be applied to study the transmission line characteristic of 
cylindrical strip lines. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
February 26, 1996 
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Node p- 1 


Figure 3. Geometry of discretized boundary curve C. 
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Figure 5. Bistatic scattering width of conducting circular cylinder excited by TE-polarized plane wave, 
a = 1 .OA, and (j) in = 1 80°. 
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Figure 6. Bistatic scattering width of coated conducting circular cylinder excited by TM-polarized 
plane wave, a = 1.0A,; b = 1.5 A,; e r = 2.0; jx r = 1.0; and c{) in = 180°. 
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Figure 7. Bi static scattering width of coated conducting circular cylinder excited by TE-polarized plane 
wave, a = 1.0A,; b = 1.5 A,; £ r = 2.0; p. r = 1.0; and (j) in = 180°. 
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Figure 8. Bistatic scattering width of coated conducting circular cylinder excited by TE-polarized plane 
wave, a = 1.0A,; b = 1.5A,; e r = 2.0; ji r = 2.0; and <|) in = 180°. 
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(a) Artificial circular boundary with radius 1 .2Xq. 

Figure 9. Isosceles triangular metallic cylinder with a = 1.85A^ and b = 1.404?Lq. FEM method is used 
inside circular boundary and BEM method is used outside circular boundary. 


19 



Y 



(b) Artificial elliptical boundary with major axis 1.3A^ and minor axis 1.1 

Figure 9. Continued. 
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(c) Artificial conformal boundary with blended comers 3 to 4 cells away from triangle. 

Figure 9. Continued. 
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(d) Artificial conformal boundary with sharp comers 1 to 3 cells away from triangle. 

Figure 9. Concluded. 
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Figure 10. Bistatic scattering cross section of isosceles triangular cylinder shown in figure 9 and 
excited by TM-polarized plane wave with angle of incidence (j)^ = 180°. Computed boundaries as 
shown in figure 9. 
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(a) Geometry and triangular mesh. 

Figure 11. Conducting strip (infinite along Z-axis) with artificial elliptical boundary with major axis 
1.3Xo and minor axis 0.3 Xq- 
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Figure 1 1 . Continued. 
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(a) Geometry and triangular mesh. 

Figure 12. Microstrip transmission line with artificial elliptical boundary with major axis 0.6 Xq and 
minor axis 0.3X 0 - 
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(b) Monostatic scattering width of microstrip transmission line excited by TM-polarized plane wave. 
W x = 0 . 9 X 0 ; W 2 = 0.1 5X 0 ; t x = 0.02Xq; t 2 = 0.1 Xq; and t 3 = O.OSXq. 


Figure 12. Concluded. 








(b) Monostatic scattering width of conducting cylinder of von Karman shape excited by TM-polarized 
plane wave. L = 2.0?io and D = L! 2. 
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